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Abstract 

A combination of small-cluster exact-diagonalization calculations and a 
well-controlled approximative method is used to examine the ground-state 
phase diagrams of the spin-one-half Falicov-Kimball model extended by the 
spin-dependent on-site interaction (J) between localized spins and itinerant 
(d) electrons, as well as by external magnetic field (h). Both the magnetic 
ordering and metamagnetic transitions are analysed as functions of h and the 
number of itinerant electrons (iV^) at selected J. Various magnetic superstruc- 
tures including axial and diagonal spin stripes are observed for nonzero values 
of J and h. Moreover, it is shown that increasing h strongly stabilizes the fully 
and partially polarized states, while the non polarized state is reduced. 

PACS numbers. : 75. 10. Lp, 71.27.+a, 71.28.+d 
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1 Introduction 



In the past decade, a considerable amount of effort has been devoted to understand 
the underlying physics that leads to a charge ordering in strongly correlated electron 
systems. The motivation was clearly due to the observation of a such ordering in 
doped nickelate [T] and cuprate |2j materials, some of which constitute materials that 
exhibit high-temperature superconductivity. One of the simplest models suitable to 
describe charge ordered phases in interacting electron systems is the Falicov-Kimball 
model (FKM) Indeed, it was shown that already the simplest version of this 
model (the spinless FKM) exhibits an extremely rich spectrum of charge ordered so- 
lutions, including various types of periodic, phase-separated and striped phases [HE]. 
However, the spinless version of the FKM, although non-trivial, is not able to ac- 
count for all aspects of real experiments. For example, many experiments show that 
a charge superstructure is accompanied by a magnetic superstructure PQI2]- m order 
to describe both types of ordering in the unified picture, it was proposed jB] a sim- 
ple model based on a generalization of the spin-one-half FKM with an anisotropic, 
spin-dependent local interaction that couples the localized and itinerant subsystems. 
The first systematic studies [UJ [7j clearly demonstrated strong cooperative effects of 
spin-dependent interaction on the ground states of the model and revealed various 
types of homogeneous as well as inhomogeneous ordering. It is naturally to expect, 
that the formation of new structures is accompanied by changes in magnetization, 
and therefore we would like to study the ground-state properties of the model in 
the external magnetic field. The motivation for our study was also the work by 
Lemahski jS], where similar studies were performed using the restricted phase dia- 
gram method for 12 periodic configurations of the spins. Here we study the model 
exactly (over the full set of spin configurations) for clusters up to 36 lattice sites and 
approximately for larger clusters. 
Our starting Hamiltonian has the form 



H = Y,tijdid ja + J^d^dn -d+d^S* 



ija i 
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where S? is a spin projection in z direction with values S? = ±1, and df a , di a are the 
creation and annihilation operators of the itinerant electrons in the (i-band Wannier 
state at site i. 

The first term of (1) is the kinetic energy corresponding to quantum-mechanical 
hopping of the itinerant <i-electrons between sites i and j. These intersite hopping 
transitions are described by the matrix elements t^, which are —t if i and j are the 
nearest neighbors and zero otherwise (in the following all parameters are measured 
in units of t). The second term is anisotropic, spin-dependent local interaction of the 
Ising type between the <i-electrons and spins. The last two terms represent energies 
of (i-electrons and spins in an external magnetic field. 

Using the fact, that S* takes only two values S* = ±1 the Hamiltonian (1) can 
be rewritten as 

if = £^+d i(T -/i££7, (2) 

ijcr i 

where = tij + a^JS* — h)5ij. Thus for a given spin distribution s = {S(, S%, Sf}, 
the Hamiltonian (2) is the second-quantized version of the single particle Hamiltonian 
t a (s), so the investigation of the model (2) is reduced to the investigation of the 
spectrum of t u for different configurations of spins. This can be done exactly over 
the full set of spin configurations or approximately over a reduced set. In this paper 
we present a combination of both methods. For clusters up to L = 36 lattice sites 
we use small-cluster exact-diagonalization calculations and for larger clusters we use 
the well-controlled numerical method described in detail in the papers [5]. 

2 Results and discussion 

In the present work we study the one (D = 1) and two (D = 2) dimensional analogue 
of the model for the spin interaction value J = 0.5 and for a wide range of magnetic 
field h (from h = to h = 0.5 with step 0.01). To reveal the finite-size effects 
numerical calculations were done on different clusters. We have found, that the 
main features of the phase diagrams in the N^ — h plane (for D = 1 as well as D = 2) 
hold on all examined lattices and thus can be used satisfactorily to represent the 
behaviour of macroscopic systems. Here we present the one dimensional results for 
L = 24, that represent the typical one-dimensional behaviour, while for D = 2 the 
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numerical results are presented for L = 36 (the largest cluster that we were able to 
consider exactly). 

Let us start a discussion of the phase diagram with a description of configuration 
types that form its basic structure in one dimension. Manifold phases entered into 
the phase diagrams are classified according to S z = J2i — Sfj and 5J = Na\ — iV^: 
the fully polarized (FP) phase characterized by \S Z \ = L, \S%\ = Nd, the partially 
polarized (PP) phases characterized by < \S Z \ < L, < |*S^| < Nd and the non 
polarized (NP) phases characterized by \S Z \ = 0, \S%\ = 0. Comparing these phase 
diagrams with ones obtained for h = [Jj, one can see rather different behaviours. 
Indeed, while for h = the magnetic phase diagrams coincide practically over the 
whole range of model parameters, the non-zero h destroys obviously this coincidence. 
Now the phase diagrams coincide only at low rf-electron concentrations where the 
ground state is the FP phase for both the ci-electron as well as spin subsystems. But 
with increasing Nd the spin and ci-electron systems are developed differently. For 
higher values of Nd the spin subsystem still prefers the FP state, contrary to the ci- 
electron subsystem where the PP state is stabilized. A different behaviour is observed 
also for h — > 0, where the spin subsystem prefers PP state, while the ci-electron 
subsystem prefers the NP state. Besides these differences the phase diagrams of spin 
and ci-electron subsystems exhibit one same characteristic and namely that very small 
changes of h can produce large cooperative changes, indicating possible metamagnetic 
transitions. Before discussing these transitions let us discuss the complete set of 
ground state configurations, that enter to the spin phase diagram. The complete set 
of ground-state configurations consists of only a few types, including 7 different NP 
phases, 10 different PP phases and of course the FP phase. Between them one can 
find different types of periodic or perturbed periodic and non-periodic configurations, 
but the most interesting examples represent NP phases. In this case all configurations 
(with the exception of = 22) have a central point symmetry. 

Let us now turn on our attention to the problem of metamagnetic transitions. To 
reveal the structure of magnetization curves we have performed exhaustive studies 
of the model for two different ci-electron concentrations rid — 1/2 and rid = 3/4. In 
Fig. 2 we summarize numerical results for both electron as well as spin subsystems. 
In the first case (rid = 1/2) the spin magnetization curve consists of two significant 
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stairs (see Fig. 2a) rrif = and rrif = 1 (or = and rrid ~ 0.5 for (^-electrons, 
see Fig. 2b). It is very interesting, that the relative small changes of h lead to the 

significant spin reorientation from the periodic arrangement {TTIITT IITTI4} 

(pif — 0) to the FP state (rrif = 1). The similar situation holds also for rid = 3/4, 
where the transition from the NP phase to the FP phase realizes through the PP 
states (Fig. 2c). In the spin magnetization curve we have observed four significant 
stairs (rrif = 0, 1/4, 1/3 and 1) corresponding to four different configuration types 

(w = [mumw^i/4 = [mmuws^i/s = [mw 3 and Wx = o ... o, 

where the lower index denotes the number of repetitions of the block [...]). The same 
structure is observed also on the <i-electron magnetization curve as illustrates Fig. 2d. 

Let us now briefly discuss the case D = 2. In Fig. 3 we present the phase diagrams 
of the model in the Nj, — h plane obtained by exact-diagonalization calculations for 
L = 36 and J = 0.5. Comparing these phase diagrams with their one-dimensional 
counterparts one can find obvious similarities. Again the phase diagrams coincide 
only at very low <i-electron concentrations where the ground state is the FP phase, 
while for increasing Nd the spin (d-electron) subsystem prefers the FP (PP) state. 
The NP ground states for non-zero h are observed only on isolated lines (points) at 
Nd'- Nd = 18, 26 for spins and Nd= 2, 10, 18, 26 and 36 for ci-electrons. 
The two-dimensional results are of particular importance since they could shed light 
on the mechanism of two-dimensional charge and magnetic ordering in doped nick- 
elate jT and cuprate [2] materials. Here we concern our attention on a descrip- 
tion of basic types of magnetic ordering that exhibits the spin-one-half FKM with 
spin-dependent interaction between ^-electrons and spins in the magnetic field h for 
D = 2. Analysing obtained numerical data, it was found that the set of ground-state 
configurations consists of only a few basic configuration types, listed in Fig. 4. 
With increasing Nj one can see a general tendency of the system to change the 
ground state from the phase separation (observed for small iV^) through special in- 
homogeneous axial stripes (intermediate N^) to the regular "chessboard" structure 
(around Nj = L). Moreover, the unusual central point symmetry observed for NP 
phases in D = 1 case persists also for D = 2 (for all NP states with the exception of 
Nd = 10 and = 22), of course now as a mirror symmetry. Again, we have stud- 
ied the magnetizations for spin as well as ci-electron subsystems on clusters up to 
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L = 100 lattice sites. Unfortunately, we were not able to obtain the definite results, 
due to finite-site effects. 

In summary, a combination of small-cluster exact diagonalization calculations 
and a well-controlled approximative method was used to study the ground-state 
phase diagrams of the generalized spin-one-half FKM with spin-dependent on-site 
interaction J = 0.5 extended by magnetic field for one and two dimensional cases. 
For both cases it was found that the magnetic field stabilizes the FP (for spins) or 
PP/FP states (for ^-electrons) , while the NP states persist only for h — > 0. Although, 
the non-zero h destroys the coincidence between the spin and <i-electron subsystems 
(observed for h — 0), the one common feature is evident, and namely, that the very 
small changes of h can produce large cooperative changes in the ground state of 
the model. This observation was supported also by studies of magnetization curves, 
where we have observed strong changes from the periodic NP/PP state to the FP 
state as well as with configuration type analysis, where we have identified phase 
separated, striped, periodic and non-periodic spin distributions. 
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Figure 1: Spin phase diagram of the model for J = 0.5 in the Nd — h plane for L = 24 
and D — 1. Inset: Magnetic phase diagram of the ci-electron subsystem. Different 
symbols represent different magnetic phases: FP (.), PP (o) and NP (+) phases. 
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Figure 2: a) The total magnetization of the spins for J = 0.5, rid = 1/2 and different 
clusters (L = 24,60 and 120) in D — 1. b) Magnetization of electrons for the same 
parameters, c) The total magnetization of the spins for J = 0.5, = 3/4 and 
different clusters (L = 24,48 and 120) in D — 1. d) Magnetization of electrons for 
the same parameters. 
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Figure 3: Spin phase diagram of the model for J = 0.5 in N d — h plane for L = 36 and 
D = 2. All ground-state configurations (with the exception of FP (.)) are displayed 
in Fig. 4. Inset: Magnetic phase diagram of the <i-electron subsystem. Different 
symbols represent different magnetic phases: FP (.), PP (o) and NP (+) phases. 
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Figure 4: Arrangements of localized spins that enter into the phase diagram displayed 
in Fig. 3. First 11 phases correspond with NP phases. All other phases are PP. To 
visualize the spin distributions we use (•) for the up spin orientation and (.) for the 
down spin orientation. 
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